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The polarization process when polarizers act on an optical field is studied. We give examples 
(3JT)' for two kinds of polarizers. The first kind presents an anisotropic absorption - as in a polaroid 

^ , film - and the second one is based on total reflection at the interface with a birefringent medium. 

■ Using the Stokes vector representation, we determine explicitly the trajectories of the wave light 
polarization during the polarization process. We find that such trajectories are not always geodesies 

ON , of the Poincare sphere as it is usually thought. Using the analogy between light polarization and 

special relativity, we find that the action of successive polarizers on the light wave polarization is 
^ ' equivalent to the action of a single resulting polarizer followed by a rotation achieved for example by 

^ , a device with optical activity. We find a composition law for polarizers similar to the composition 

* J~i . law for noncollinear velocities in special relativity. We define an angle equivalent to the relativistic 
f) l' Wigner angle which can be used to quantify the quality of two composed polarizers. 
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^ ■ I. INTRODUCTION 

It is now usual to describe the polarization state of an electromagnetic field by using the notion of Stokes vector 
and that of the Poincare sphere. This elegant representation of the polarization state is very useful in optics, e^g 
ly-j ' geometrical phases such as the Pancharatnam phase can be efficiently interpreted on the Poincare sphere jl| [|. 

, Actions of polarizing devices are described in such a space with simple mathematical operations. As an example, 
C*") ■ variations of the polarization state when light passes through an optical active medium is expressed by a simple 
rotation around the axis connecting the two poles of the sphere. In the same way, the effect of a birefringent plate 
qq ■ corresponds to a rotation around a vector lying in the equatorial plane by an angle determined by the optical path 
delay between the ordinary and the extraordinary axis of the plate. However, it is curiously interesting to note 
that the evolution of the polarization due to the action of a polarizer has never been studied in details. On the 
Poincare sphere, such an evolution corresponds to a trajectory intuitively supposed to be a geodesic connecting the 
polarization states of the field before and after the polarizer (see e.g. Q or ||). As we will see, the study of the 
actions of two different kinds of polarizers (the first one using anisotropic absorbing medium and the second one 
using total reflection of one component of the electromagnetic field at the interface with an anisotropic media) will 
show that such an assumption is not always true and that the trajectories of the polarization state of the field can 
be more complex and do not necessarily correspond to Poincare sphere geodesies. At the beginning of the article, we 
introduce the Jones representation of a polarizer. The action of a polarizer on the polarization state is then studied by 
introducing rotation operators and Lorentz boots operators borrowed from special relativity. Since the SL(2, C) group 
is homomorphic to the Lorentz group (see e.g. Q), a formal equivalence between the special relativity space and 
the polarization space exists S, [|. This formal equivalence can then be used to determine what are the counterpart 
in the polarization space of, for example, the composition law of velocities, of the aberration phenomenon or of the 
Wigner angle. 
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The paper is organized as follows. In section II we define the mathematical formalism used to describe the effect 
of polarizers on an optical field polarization. The states of polarization of the field are represented on the Poincare 
sphere and the action of the optical devices is defined by their associated Jones matrix. It appears that the effect of 
polarizers on the polarization are essentially defined by a complex number A which contains informations about the 
propagation and the absorption of light in the polarizer. In section III we analyze the evolution of the Stokes vector 
s when light goes through different kinds of polarizers. Examples are studied and the trajectories of given optical 
states during the polarization process are analyzed. We point out that the projection of such trajectories on the 
Poincare sphere are no longer a geodesic of that sphere as it could have been classically expected. In section IV we 
focus on the evolution of the optical field intensity and we derive a natural generalized Malus law. Then we examine 
the action of two successive polarizers and we show that they are equivalent to one polarizer followed by a rotation 
exactly as in special relativity the composition of two Lorentz boosts is a Lorentz boost followed by a rotation. In 
section V we find the general composition law which permits to determine characteristics of the equivalent polarizer 
and the angle of the rotation. Then we discuss the interpretation of such a composition law. In section VI we give 
the counterpart of the Wigner angle in special relativity for two composed polarizers and we relate such an angle to 
the quality of the polarizers. 



II. MATHEMATICAL DESCRIPTION 



A. Stokes parameters and the Poincare sphere 



The polarization state of light will be described with the use of the Poincare sphere. Our goal is to describe the 
evolution of the polarization state during its propagation inside a polarizer. Without loss of generality we consider 
that the optical field propagates along the z direction. This field can then be described by the Jones vector 
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where E x and E y are the two complex components of the electric field. Instead of the basis {\x}, \y}} we prefer to 
work with the basis {\R), \L)} related to right- and left-handed circularly polarized states 
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In this basis the Jones vector \ip) ([!]) is transformed into 



i 

72 



Ev 



iEy 

i E„ 



y / {\R),\m 



(2) 



(3) 



The intensity I and the polarization of the wave light can be characterized by the Stokes parameters p 
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If we construct a three dimensional vector ~s (called hereafter Stokes vector) with the components s , s 2 , and s 3 , it 



is easy to check that its norm s 
using spherical coordinates as 



/. Thus, for a given intensity, the Stokes vector can be parametrized 



s = s | sin 8 cos (jTe 1 + sin 9 sin d>~~e 2 



cos ft e 



(5) 



where 6 € [0,tt], € [0,2tt[ and {e , e ,~e 3 } is an R 3 orthonormal basis. The Stokes vector ~s (|J) determines 
a point located on the Poincare sphere S 2 of radius s. Since the norm s gives the intensity of the wave light, the 
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Right handed circularly 




polarized states 

FIG. 1: Representation of the Stokes vector s and the Poincare sphere S 2 

direction of the vector ~s, i.e. the angles 9 and <fr, characterizes the polarization. Thus, the polarization state of the 
light wave corresponds to a unique point on the Poincare sphere S 2 . 

In order to follow the evolution of the light wave polarization, i.e. the evolution of the Stokes vector ~s , we prefer 
to work with the projector \ip) instead of \rp) itself. Indeed, this projector can be easily expressed as a function of 
~s since 

mm = K *° x J i - i f ) = \ ^ + ~* • s ^- (6) 

Here cr° is the 2x2 identity matrix and ~a is a vector the components of which are the usual Pauli matrices 

M°»)< - 3 =0-°)- < 7 » 

Usually |Tc[] , one choose as the north (south) pole of the Poincare sphere the circular right-handed (left-handed) 
polarization state. With this definition, the angles 9 and 4> (||) are directly related with the ellipticity angle \ an d 
the azimuth angle \t since 9 = ir/2 — 2\ and <f> = 2\f respectively [0. The circular right(left)-handed polarization 
state corresponds to 9 = (9 = ir). Right-handed (left-handed) clliptically polarized states correspond to s 3 > 
(s 3 < 0), or cquivalcntly x > (x < 0). Linear polarization states are represented by Stokes vectors lying in the 
equatorial plane of the Poincare Sphere since 9 = ir/2 (x = 0). Each linear polarization states is determined by a 
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given angle </>, or equivalently by a given azimuth angle '5. For a linear polarized wave, the angle $ is defined as 
the angle between the vibration plane and a laboratory axis e x (reference axis) orthogonal to the direction of the 
propagation. The difference of polarization for two linear polarized light waves 1 and 2 can be measured, in the 
laboratory, as the difference A\l/ = — $2, or equivalently, on the Poincare sphere, as A(j> = <pi — <p2 ~ 2A\f f . Thus 
an angle measured in the laboratory is half of the corresponding angle in the Poincare sphere representation. For 
example, two orthogonal polarizations are characterized by A "J = 7r/2 in the laboratory frame and by A<p = 7r in the 
Poincare sphere representation. Although the rest of the paper does not consider any particular choice of basis nor 
any particular choice of polarization, it is useful to keep in mind these last remarks in order to be able to make at any 
time the parallel with known conventions in optics. So, in a general manner, two states with orthogonal polarizations 
are represented by opposite Stokes vectors, and correspond to projectors p-g and p_-j. 



B. Jones matrices for polarizing devices 

1. Jones matrix for birefringent systems 

The action of any optical system on the light wave state ([j]) is defined by a 2 x 2 Jones matrix. For example, 
the Jones matrix associated with a birefringent system and acting on l^) can be written as 

e# ° ° ^ - e ^+«/ 2 ( 6#/2 ° , ^| (8) 
e'*« )~ e { e"^/ 2 ) ■ ( *> 

Here 4> and 4>e are the ordinary and extraordinary phases associated with the propagation of the optical field polarized 
along each optical axis. We have defined the difference between these two phases as (j) = <f> — 4> e . For ordinary and 
extraordinary axis rotated by an angle 9 in the (x, y) laboratory frame, the Jones matrix (H) becomes 



l( o +0 c )/2 f cos6 sinO \ f e^' 2 



B(6, <j> , M = e«W ^ — j ^ Q J H2 ) [ Z e ~oZ ) • ( Q ) 

The corresponding Jones matrix for a birefringent optical device acting now on the Jones vector \ip) written in the 
{\R),\L}} basis (§ is then 

B{6, 4> , e ) = \( K \ -\\ B{6, <p o , 4> e ) (] _] V (10) 
After performing the matrix multiplications, it is possible to recast (fol) in terms of Pauli matrices 



B(6, cj> , e ) = e ^° + ^' 2 [cos p + isin^.Vj (11) 

where ~p = cos 201? 1 + sin20~e >2 . Using now the mathematical relation e za q ' a = cos a + isinaTq ■ a valid for any 
angle a and any unitary vector ~q ', we can rewrite ([ll]) in a more compact manner 

B(6,4> o ,(t))=e i 't >0 e- i %e i % ir ^. (12) 



As we will see in section III A, the action of the matrix B(e, (f> , 4> e ) on the polarization state \ip) can be viewed in 
the Poincare sphere framework as the rotation of the Stokes vector ~s around the vector ~p by an angle </>. As no 
dissipation phenomenon occurs, the norm s° of the vector s is constant. 



2. Jones matrix for polarizers 
The Jones matrix associated with a perfect polarizer acting on can be written as 
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In that example only the x component of the field is conserved, the orthogonal component is not transmitted. This 
oversimplified picture corresponds to the limit case where the attenuation factors in the x and y directions differ from 
each other by several order of magnitude. The realistic case correspond then to the matrix 



e -(7 1+ , 2 )/2( - U ) (W) 



" 71 \ _ ( 7l+72 )/ 2 ( e 7/2 
e-^ ) \ e" 7 / 2 

where e -71 and e~ 72 are the attenuation factors in the x and y directions respectively. We have defined the difference 
between the two attenuation terms as 7 = 72 — 71 ■ Again, if the polarizer axis is rotated about an angle , in the 
(x, y) laboratory frame, the Jones matrix ( fl4| ) becomes 

P(,, 7l , 72 ) = e-^y, ( cos* sin, W e 7 \ / cos, -sin, \ _ 



sin, cos, y I e 7 / I sin, cos 

Following the same steps as in the previous section for birefringent systems, the Jones matrix for a polarizer acting 
on \ip) in the {\R), \L)} basis (§) is 

P(,, 7l ,7) = e" 7l e-*e*P > - 1?> (16) 

where the vector ~p encodes again the information about the direction of the polarizer in the laboratory frame. In 
the Poincare sphere framework, as the attenuation difference 7 increases, the action of the operator (um on can 



be seen as the progressive collapse of the Stokes vector s on the polarizer vector p direction (see section [II B). 

We now use the derived expressions (1^) and ( |l6| ) to examine more precisely the case of a polarizer based on 
anisotropic absorption (such as a polaroid film) , and the case of a polarizer based on a total reflection at the interface 
of an anisotropic crystal. 

The case of an anisotropic absorbing polarizer can be modeled by a system where the phases <fi Q and <j) e in orthogonal 
directions are complex, 

4>o = (k + ia a )z, , . 

4> e = (k + ia e )z. 

The parameter k is the wavevector of the optical field in the film and a and a e are the two absorption coefficients 
(for the amplitude) in two orthogonal directions of polarization. The system acts as a linear polarizer if a ^ a e . In 
equations (|l7|), the parameter z is the propagation coordinate, i.e. z = at the entrance of the polarizer and, e.g. 
z = L at its end. Hence, using ([T2|), the operator acting on \ip) associated to such a polarizer is 

Pi(,,</» , 7 ) = e tkz e- a - z e~i (18) 

where 7 = (a e — a a )z. Up to a global phase term and a global attenuation term which do not affect the light wave 
polarization, we retrieved, as expected, the expression of the Jones matrix for a polarizer (|Tc|). 

The other case corresponds to a polarizer based on total reflection. We consider a polarized light hitting the face 
of an anisotropic crystal in such a manner that only one component of the polarized light is transmitted, the other 
component being evanescent after the interface. This kind of polarizer can be modeled by an optical device with the 
phase <j) and <j> e defined as 

4>o = kz, 

(fi e = iflZ 

where k and /i are the z real components of the wavevector of the linear polarization which are transmitted and 
totally reflected respectively. Here, the parameter z is the propagation coordinate along the transmitted light wave 
component (z = at the interface). Again, using (p~2|) . the Jones operator acting on the wave light and associated 
to such polarizer is 

P 2 (,,5, 7 ) = e ^ e -^ e^^- V (20) 

where 7 = fiz and 8 = kz. 

Comparing the just derived expressions ((lj), (pit), ( |l8| ) and ( po|) of Jones operators (Jones matrices) for polarizing 
devices, we remark that all of them are written in the following manner 

P(p,<i> p ,A) =e i ^e~^e^ (21) 
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where A is a complex differential attenuation term, <f) p an overall phase/attenuation term, and the polarization 
vector associated with the axis of the polarizer. In the following we will drop out the overall factor e l ^ p from the 
expression of the Jones operator (|2l|) since it docs not affect the evolution of the wave light polarization. We have 
just to keep in mind that in the case of linear polarizer such as those defined in (H) and ©, the overall factor 
contains an attenuation part which globally reduces the intensity of the wave light. From now on, it will be useful 
to split A in its real and imaginary parts, A = 7 + i 6. Note that A, and consequently 7 and 8, are implicit linear 
functions of the z propagation coordinate. 

Also, although equation (|2l| ) has been derived from the analysis of linear polarizers (vector ~p lying in the equatorial 
plane), our model is completely general. It is possible to use a more general rotation than ( |l5| ) bringing the polarizer 
characteristic vector ~p outside the equatorial plane. For example, if the vector ~p points toward the Poincare sphere 
North (or South) pole, our model pl|) describes circular polarizers or, if A = iS, media with optical activity. 



III. LIGHT WAVE POLARIZATION VIEWED AS A STOKES 4- VECTOR TRANSFORMATION 

According to the previous section, the action of a polarizing device on a light wave can be characterized by the 
following operator 

PHnA z ) = e-^ +l5) e^ +l5) ^^ , (22) 

where 7(2) and 5 (z) are real positive functions of the transformation parameter z £ R + . Usually, z is the length 
penetration of light into the device and the functions 7(2) and S(z) encode the absorption and the propagation of 
light inside the device respectively. We restrict our work to homogeneous media for which these functions are linear 
in z, i.e. 7(2) = az and S(z) = (3z with a, f3 £ R. The vector ~p is a normalized polarization vector (|| ~p\\ = p = 1) 
associated with the polarizing device. The vector p determines a point on the Poincare sphere which corresponds to 
a pure state of polarization p . 

Using the projector = we are able to define the Stokes 4-vector s = (s , ~s) with the following components 



s" = TrKp t ), M e {0,1,2,3}. (23) 

It is interesting to note that, introducing the Minkowskian metric 77 M „ = diag{— 1, 1,1,1}, the Stokes 4-vector can be 
considered as a light-like 4-vector since rj^vS^s" = 0. The spatial components of the Stokes 4-vcctor, {s z }i=i...3, are 
the components of the three dimensional Stokes vector ~s defined in the previous section. The time component of 
the Stokes 4-vcctor s° is the norm of the Stokes vector s . Consequently this last component also corresponds to the 
light wave intensity s° — (ip\ip) = I. 

So, according to the operator (p2[), a polarization state />?(o) is continuously transformed into another polarization 
state PY(z) by the following operation 

^(*)=-P^,7,ff(*)^(0) i ^ l7l4 (*). ( 24 ) 
Now, using ( p3] ) and (^) we are able to define the evolution of the Stokes 4-vector s(z) by 

8»{z) = Tr( CT ^ (z) ) 



Tr (p^ 7 j (z)<T"F^ 7 j (z)^ (0) ) , fx £ {0, 1, 2, 3} 



(25) 



where we recall that s°(z) = s(z) = \\ s (z)\\ = I(z) is the intensity of the light wave, and {s l (z)}ie{i,2.3} are the 
spatial components of the Stokes vector s (z). After some algebraic calculus, the components s^{z) of the Stokes 
4-vector (G5) can be expressed as 



s°(z) = e"T( z ) (s°(0)cosh 7 (z) + sum 7(2)^(0) • ~p ) 

s "( z ) = { ^>( z ) = e - 7 («) ( ^>(o) cos 6(z) + ~s(0) A ~p sin S(z) + (1 - cos 6(z)) {~p ■ 1>(0)) ~P ^ 

+ s°(0)"p sinh7(z) + (cosh7(z) - 1) (~p ■ ~s*(0))~p ) 

For the sake of clarity we now consider the two particular cases 7 = , 5 7^ and 7 7^ , 5 = before considering the 
more general case of polarizers with 7 7^ and 6^0. 



FIG. 2: In the case of a non absorbing polarizing device, the Stokes vector s is rotated around the polarization vector p. 



A. Non absorbing polarizing device : 7 = 0, 5 7^ 



The unitary case 7 = 
activity. In such a case, 



0, S =/= 0, \/z G R corresponds to non absorbing birefringent devices or to media with optical 
the transformation (Eq) of the Stokes 4-vector s becomes 



,0 



(z) = S °(0) = 1 



~s{z) = ~s (0) cos S(z) 

+ ~s(0) A ~p sin<5(z) 

+ (1 - cos5(z)) ("p • ~s'(0))~p. 



(27) 



The first equation in (27) verifies that s(z) = || s (x)|| = s (0) = 1 is invariant under the action of the non absorbing 
polarizing device. So, the light wave intensity / is conserved and for convenience we have set 1=1. The extremity 
of the vector ~~s(z) still lies on the Poincare unit sphere after the transformation. The second equation in ( |27| ) is the 
Rodrigues Formula for the rotation. This equation implies that the Stokes vector ~s(0) has been rotated around the 
vector ~p by an angle S(z) (Fig. ||). Without any loss of generality, we can always choose an K 3 orthonormal basis 



{e}, 



ie{l,2,3} 



such as 



P 



Then the Stokes 4-vector transformation (J27|) can be written in the following simple 
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form 

/ 1 

s ( z ) = n n 'r/ \ . ! n I s(0) 
w 10 coso(z) smo(z) 1 

\0 — sin<5(z) cosS(z) 

where we recognize the rotation matrix around the "e 1 axis by an angle 5{z). 

The transformation ( |27| ) of the light wave polarization presented in Fig. |^ is typical of birefringent devices. For 
example a A/4 plate can be used to transform a linear polarized light wave into a circular polarized one, i.e. to 
transform a Stokes vector lying in the equatorial plane into a Stokes vector pointing toward a pole of the Poincare 
sphere. 

B. Absorbing polarizer : 7 7^ 0, 8 = 

We consider now the non-unitary case 7 7^ 0, S = 0, Vz £ K+. In such a case, the transformation (|2^) of the Stokes 
4-vector s becomes 

f s°(z) = (s°(z) cosh 7 (z) + sinli7(z) V(0) • ~p) 

= { (29) 
[ = e-TW (-?(0) + s°(0) "p sinh7(z) + (cosh 7 (z) - 1) ( ~p ■ ~f(Q)) ~p) 

Naturally, the s° component giving the intensity of the light wave (first equation in (pjj)) is a decreasing function of 
the parameter 7. As 7 — ► 00, the Stokes vector is brought along the direction of the p vector since 

~s(z) 7 oc ~p . (30) 

For a finite value of 7, the Stokes vector is not completely brought along the direction of the vector p , the polarization 
of the light wave is then only partial. From now on, devices which transform the light wave polarization according to 
( p9f ) will be designated as non perfect polarizers (finite 7) or perfect polarizer (7 — > 00). 

1. Four dimensional representation 
Again, without any loss of generality, we can choose ~p = it 1 . Then the transformation reads 



s(z) = e 



-7(«) 



/ cosh 7(2:) sinli7(z) 

sinh7(z) cosh 7(2) 

10 

V 1 



s(0). (31) 



We recognize here a Lorentz boost in the ~p — ~e 1 direction, the absorption term 7 is then similar to the rapidity 
$ = argtanhw/c in special relativity. So, up to an attenuation factor e -7 , the transformation ( |29| ) is a Lorentz boost 
in the direction of the polarizer vector p . The Lorentz boost mixes only the intensity component s°(0) and the 
component along the polarizer vector, i.e. the projection s (0).p . The other components are left unchanged by the 
boost. Associated with the global attenuation e~ 7 ^ z \ the boost attenuates the components of V which arc orthogonal 
to the polarizer vector ~p: ~s is progressively brought in the direction of the polarization vector ~p (Fig. ||). 
Note that this case exactly corresponds to using a polaroid film with 7 = (a e — a a ) z. 

2. The Poincare sphere representation 

The evolution of the light wave polarization can also be followed on the surface of the Poincare sphere by studying 
the evolution of the normalized Stokes vector s (z) = ~~s(z)/ ||~?(z)|| 
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FIG. 3: (Color online) Transformation ([n]) of the spatial components s(z) of the Stokes 4-vector s(z). The polarization vector 
is p = e 1 . The incoming light wave has an initial polarization s(0) = l/\/3 (e* 1 + e 2 + e 3 ) represented by the thick red (light 
gray) dashed vector. The Stokes vector for 7 — > 00 is alignated along the polarization vector p = e 1 and is represented by the 
thick blue (dark gray) vector. The other vectors are drawn to show the intermediate positions of the Stokes vector s as the 
parameter 7 (or equivalently the parameter z) increases from to 00. 



a"(z) = < 



liz) 



s°(0) = s°(0) = 1 

~s(0) + ~p sinhj(z) + (cosli7(z) — 1) ( p ■ 1^(0)) p 
s°(0) cosh7(z) + sinh7(z)lf(0) • ~~p 



(32) 
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FIG. 4: (Color online) Representation of typical evolutions of the Stokes vector s(z) and of the normalized Stokes vector s(z) 
according to the transformations ( p9[ ) and ( p^ ) respectively. 

We remark easily that throughout the continuous transformation the Stokes vector ~s(z) is always orthogonal to the 
unit vector 

g= , m 1 m ?At(Q) (33) 
s u (0) sma(0) 

where a(z) = ("p*, ~s*(z)). Indeed, for any parameter z, we can write 

~s(0) A ~s*(2) oc ~q. (34) 

The set of the Stokes vectors ~~s(z) with z G K + defines a plane orthogonal to the vector ~~q and containing the center 
of the Poincare sphere. This implies that the normalized Stokes vector s (z) describes an arc of a great circle of the 
Poincare sphere S 2 (Fig. ||). The projection of the Stokes vector ~s(z) describes a geodesic of the iS 2 unit Poincare 
sphere. Therefore, the normalized Stokes vector s (z) is rotated around the axis q by an angle d(z) = a(0) — a(z) 
(see Fig. ^). Using the Rodrigues formula, we can rewrite the transformation ([32]) as a simple rotation 

r s°(z) = s°(o) = i 

s»(z) = \ (35) 

[ JL( z ) = ~?(0) cosd(z) + ~~s(0) A ~q sinS(z) 

with 



sin 8{z ) = sina(0) cosa(2;) — sina(z) cosa(0) 



(36) 



11 



and 

cosi5(z) = cosa(O) cosa(z) + sina(z) sina(O) (37) 

where 

tanh 7 (z) + cosa(0) 

cosa (z) = - — — — 38) 

v ; 1 + tanh7(z)cosa(0) 



and 



/ \ sina(O) .„„. 

sina z = — - — — . i ' 77T7. (39) 

cosh 7(2:) + smn7(z) cos a(Q) 



C. General polarizer : 7 7^ 0, S 7^ 

Choosing again ~p = ~e , the transformation (|26|) for a general polarizer can be written in the following matrix 



form 



s(z) 



„-7« 



/ cosh 7(2) sinh7(z) 

sinh7(z) cosh.7(z) 

V 



/ 1 

10 

cos<5(z) sm<5( 

\0 — sin<5(z) cos<5( 




(40) 



Equation ([40|) shows clearly that the action of a polarizer on a light wave can be seen as three successive operations: 
a rotation of the Stokes vector 1^(0) by an angle S(z) around the polarizer vector a Lorentz boost in the direction 
of the polarizer vector ~p , and the global attenuation e -7 ^. Each of these three operations evidently commutes with 
the others. Thus, the operator (123) associated with a general polarizing device can be rewritten as 



Pjt n A z ) = e ^ 7(z) e-^ (z) %, 5 (z)%, 7 (z) (41) 
where in the SX(2,C) group representation 

B 7 >Jz) = e 2 ¥?-° (42) 
is the Lorentz boost operator in the direction ~p with the rapidity 7 and 

= e iS &~** (43) 

is the rotation operator with the rotation vector ~q and the angle of rotation S. Fig. ^| gives a typical illustration of 
the Stokes vector evolution for a general polarizer (|2^). 

Note that this case corresponds to total reflection based polarizers ( p0| ) with 7 = /iz and S = —kz. 

As a conclusion of this section, we can state that during the polarization process of a wave light, the evolution of the 
Stokes vector, i.e. the evolution of the polarization, does not necessarily describe a geodesic on the Poincare sphere. 
A geodesic, i.e. here a part of a great circle of S 2 , is found only in the case of polarizing device such as a A/4 plate 
or in the case of non-unitary polarizer with 5 = such as e.g. a polaroid film. In this last case, it is the projection of 
the Stokes vector on the Poincare sphere which describe a geodesic. All these statements on the effective trajectory 
of the light wave polarization are worth to consider when geometric phase have to be calculated O] . 

IV. GENERALIZED MALUS LAW, DEGREE OF POLARIZATION 

Now, let us look at the intensity component s°(z) of the Stokes 4- vector. We define a(z) = 29(z) as the oriented 
angle (^p,~s(z)) between the polarizer vector ~p and the Stokes vector ~s{z) for the parameter z. Then, 9{z) is the 
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FIG. 5: (Color online) Typical representation of the Stokes vector s(z) evolution according to the general transformation ( |4C| ) 
or (|2(]). The two dashed lines give the direction of the initial Stokes vector s(0) and the direction of the polarization vector 
p respectively. The initial Stokes vector s(0) is represented by the thick black dashed vector. The Stokes vector for 7 — > 00 
is represented by the thick black vector oriented along the direction of the polarization vector p. The other vectors are the 
successive Stokes vectors s as the parameter 7 (or equivalently the parameter z) increases from to 00. 



angle, in the laboratory frame, between the axis of the polarizer and the ordinary axis of the light wave at position z. 
The intensity of the light wave, given by the first equation in (26), can be rewritten as 



I(z) = 1(0) (cos 2 9(0) + sin 2 6(0) e - 2 ^ z ^ 



(44) 
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a(0)=26(0) a(0)=29(0) 

FIG. 6: (Color online) (a) Generalized Malus law: representation of the ratio I(z)/I(0) as a function of a(0) = 29(0) (angle 
between the polarizer vector p and the incoming wave polarization s(0)) for different values of the absorption term f(z). (b) 
Angle a(z) between the polarizer vector p and the wave polarization s(z) as a function of a(0) for different values of the 
absorption term 7(2). 



which is a generalization of the Malus law for non perfect polarizers. Fig. [j]a shows the Malus law for different value 
of the absorption term 7(2) from to +00. Let us define I\\(z) = I-^(z) as the intensity of the light wave in the p-^> 
polarizer state and I±(z) = I--jf(z) as the intensity in the corresponding orthogonal state p--g. After some calculus 
we obtain 

I ]l (z) = I^(z) = Ti{p^ (z) p^) =7(0) cos 2 0(0) (45) 

and 

I ± (z) = I_ r {z) = Ti(p^ (z) p_^) = /(O)sin 2 0(O)e- 2 TW (46) 
Now, if the polarizer is placed in an unpolarized beam the transmitted intensity at the point z is 

(m)e(o) = (W)o<f» + fr-(*)>fl ( o) = ^ (l + e ' 27(2) ) (47) 

where the brackets (• • • }gm^ denote the average over the initial angles 9(0). From (pTj). we remark that for a non 
perfect polarizer (7 finite) more than the half of the incoming intensity is transmitted. The degree of polarization 
r](z) of the non perfect polarizer is given by 

(W) 9 (0)-(^(*)) 9 (0) 

»7(*) = 77TT\ Z77 n\ = tanh 7 (z) (48) 



and the extinction ratio £(z) by 



fW = ^n^=e-^W. (49) 



As 7(2) varies from to +00, 77(2) and £(z) vary from to 1 and from 1 to respectively. Although rj(z) and £(z) 
measure already the efficiency of the polarizer, we can always ask how far is the polarization vector ~s (z) from the 
polarization vector p associated with the polarizer. For that purpose we focus on the angle a(z) between the two 
vectors. Combining the two equations in (26), or equivalently using ( f48|) and (|38|), we obtain the following relation 

r)(z) + cosa(0) 

cosa(z) = — - — — — -. (50) 

w 1 + T)(z) cos a(0) 

from which we arc able to extract the angle a(z) as 

a(z) = arccos (tanh (7(2;) + argtanh (coscv(O)))) . (51) 
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FIG. 7: Two successive boosts with respective direction pi and p*2 can be decomposed in a boost of direction p and a rotation 
around the direction q. 



Thus, as 7(z) varies from to +oo, the angle a(z) varies from a(0) to 0, or equivalently cosa(z) varies from cosa(0) 
to 1. Fig. shows a(z) as a function of a(0) for different values of the absorption term j(z). Expressions ( |50| ) and 
( ^l| ) are dependent of the initial angle a(0). The integration of ( ^p|) over all the possible initial angles a(0), i.e. 

i r 

(cos a(z)) IQ) = — da(0) cos a(z), (52) 
K ' 71 Jo 

provides a convenient parameter to characterize the quality of the polarizer. Analytic calculation of (E^) gives us 

(00Bo(*)> a(p) = 1 ' V/ ?? 1 ( ~ ??2(Z) = tanh^. (53) 
Hence, as j(z) varies from to +oo (perfect polarizer), the parameter (cos ce{z)} a / \ varies from to 1. 



V. COMPOSITION LAW FOR POLARIZERS 



From now on, we assume the z parameter dependence of the different functions, such as e.g. S or 7, and for clarity 
we drop the z parameter from the following mathematical expressions. 
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A. Association of two polarizers 



Let us now consider the successive actions of two absorbing polarizers (S = 0, see section IIIB) with different 
polarization axis p 1 and p 2 , and different absorption terms 71 and 72. The corresponding operator of the total 
system is then 



P 2,72,0-* p i,7i,0 



p -(7l+72)/2 r> D_ 

P 2,72""! 



i,7i 



(54) 



where we use the SL(2, C) group representation ( f42| ) of the Lorentz b oost operator. We keep in mind that the case of 
the composition of two general polarizers, i.e. with 5 7^ (see section III C ), can be easily retrieved from the case (H) 
of two absorbing polarizers (8 — 0). We focus now on the two successive boosts in (K3). We know (see e.g. P, M) that 



the action of two successive noncollinear Lorentz boosts, By x 7l and -Bp> 2j72 , is equivalent to the action of a Lorentz 



boost B-^ n followed by a rotation R-f 5, i- 



P 2,72 P 1,71 



' g ,(5-° £> ,7 



with suitable parameters 7 and 8, and suitable directions ~p* and ~q . Then, equation ( pi| ) can be rewritten as 



P— > P— > 

P 2,72,0-* p i,7l,0 

or, using only polarizing device operators (p3), as 



„-(7l+72)/2 n d 

^ *- g ,0 p ,7 



— p(7-7i — n)l^ p iS/2p p 
P 2,72,0-' p 1,71,0 0-1 g ^jd- 1 p ,7,0- 



(55) 



(56) 



(57) 



In the appendix, we derive from the initial parameters 71, 72, ~p 1 and ~~p 2 the expressions of the resulting rotation 
parameters 8 and ~~a and of the resulting boost parameters 7 and ente ring the relation ( |55"| ) and consequently 
entering equations (p6| ) and ((5?]). For the rotation parameters we have (|7~ 



1 



sin a 



P2A pi 



and (A.c) 



5 = 2 arg (l + tanh y tanh y e M ) 



(58) 



(59) 



where a = (~pi,~P2) € [0,tt] is the angle between the two polarizer vectors ~p 1 and ~p - 2 - Equation ( pq) imply that 
the rotation axis ~(f is orthogonal to the plane defined by the vectors p \ and ~~p 2- Since the vectors ~pi, ~~Pi and 
~p belong to the same plane (see appendix section Q]) , it is possible to parametrize these vectors with the help of an 

arbitrary reference axis u belonging to this plane. Hence, the angles a\ = (1?, p 1), a-i = (1?, p 2), and <j) = ( u , ~p) 
characterize completely the respective vectors ~p\, ~p 2 and ~p. F igure jj] gives an illustration of the relative positions 
of these vectors and the associated angles. According to ( A. 14 ), the boost parameters are given by the following 
expression 



7 

tanh — e l 



tanh — e i 



tanh — e ia2 



1 



tanh ^ tanh ^ e l( " 2 



-Ul) 



(60) 



Indeed, the modulus and the argum ent of 
resulting polarizer (through (j> ( A. 16| ) ) . 

Let us now define the following addition law 



2 2 

) give respectively the absorption term 7 and the direction ~p of the 



a®b = 



a,b e 



Using that law, equation (r 



l + a*& 

can be rewritten in the following more elegant form 



Tie" 1 



r 2 e 4 



where we have have noted T = tanh % . Since 



P— > P—> or R— ► 

P 2,72,0-* p 1,71,0 ^ P 2,72 P 1,71 



Q ,5 p ,7 



(61) 



(62) 



(63) 



equation (|6^) , which is similar to the addition law for non collinear velocities in special relativity jt4| , is thus also the 
composition law for polarizers. 
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B. Association of TV polarizers 



Consider now the action of N successive absorbing polarizers (see section [II B). The corresponding operator is 

V(1,2,...,N) = IWV ii7 , i0 

= P_ ...P. P_ 

Pn.7jv:° P 2 >7 2 .° Ph7i>0 

/ N \ ^ N 

= expf-lEf-r^)^ fCV. (64) 

= exp (-1 Ef=i 7,j ^ - 1, /V)S(A - 1, JV)II 3 =i ^., 7 . 

= cx P (-i Ef=i 7i) ( r£iW ■ ■ ■ > J ■ ■ ■ 



where we use the following notations 



I', .. «„•.. K(i,j)B(i,j) (65) 



and 



B(i, N)B t =TZ(i-l,..., N)B(i — 1, . . . , N). (66) 



The operator B(i, . . . , N) is the boost operator resulting from the combination of the (N — i) last boost operators. The 
operator 7Z(i, . . . , N) is the rotation operator resulting from the combination of the boost operators B(i + 1, . . . , N) 
and B-, . For the sake of completeness, we define B(i) = B^ ^ and H(i) = a . In (p^), the boost operator 

B(l J ... ) JV) = B^ 7 =e^-" ? (67) 

resulting from the combination of the iV boost operators B^ (j = 1, . . . , N) is characterized by an absorption 

p j ' 7 j 

parameter 7 and a vector p directly calculated from the N — 1 successive applications of the composition law (|62j). 
The ordered product in (|64|) 

J] K(j,...,N)=R^ =6*1* (68) 

is a resulting rotation operator characterized by a rotation vector if and an angle 5 calculated from successive 
applications of Eq. ( A.9). The operator ( p4[ ) corresponding to N successive absorbing polarizers can then be rewritten 
as 

V(l,2,...,N) = e-^j=^iR B-> (69) 

q , S P il 

or using polarizing device operators ( p2|) as 

V(l,2,...,N) = e3^-^f=i^)e* 5 P_^ n . (70) 

v ' ' ' y "9,0,5 p. 7.0 v 1 

Using ( |69[ ) , the initial state p^> of an incoming light wave is transformed into the following final state 

p r =P(l,...,N)p lto P\l,...,N) (71) 

with the corresponding intensity 

I = Tr(p-) = /oe 7 "^ 1 7j (cos 2 0„ + sin 2 6> e- 27 ) . (72) 



Here 29q (sec section IV ) is the angle between the polarization s of the incoming wave and the resulting polarization 
vector p . So, the transmitted intensity of an unpolarizcd light beam is 

W e „ = T e7 " E " l7j ( 1 + e " 27 )- ( 73 ) 
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Let us define r\-g as the degree of polarization of the N polarizers using the resulting polarization state p-g as the 
state of reference. Using, 

I t = Tr(p r pf) = JoeT-E^ ^ cos 2 6 (74) 

and 

I_f = Tr( Plt p_ t ) = J e-T-^ 7 ' sin 2 O (75) 
we obtain the corresponding degree of polarization 



and the corresponding extinction ratio 



<%> 9o + ( 7 -P> eo 



tanh7 (76) 



e" 27 . (77) 



P /8o 



The expressions found in ([76]) and ^77]) arc similar to the expressions of the degree of polarization r\ ( ^8|) and of the 
extinction ratio £ ( |49| ) found for one polarizer except that here the absorption term 7 results from the action of N 
absorbing polarizers. So, the composition law (|62|), whose the N — 1 successive application give T = tanh ^ is also 
the composition law for the degree of polarization rj (|76|). 

To conclude this section let us retrieve the intensity of an initially unpolarized light beam through the succession 
of two perfect polarizers (71 — > 00, 72 — > 00). First, the modulus of the composition law (60) implies that 7 — ► 00. In 
that case the expression ( |73|) becomes 

(/)„ = ^ e 7-7i-7*. ( 78) 



Using equations ( |A.2| ) and (A. 8) in the case of two perfect polarizers (71 — > 00, 72 — -> 00, 7 — ► 00) we obtain the 
following relation 

cos 2 I = e 7 " 71 " 72 ( 79) 

where a = 2(9i2 is the angle between the polarization vectors p 1 and p 2 (in the Poincare sphere representation) of 
the two successive polarizers. Thus Eq. ([78]) can be rewritten as 

(J) 9o = |cos 2 ^ 12 (80) 

which is the well known expression for the intensity of an initially unpolarized light beam going through two perfect 
polarizers. We remark that 7 given by the composition law (|m]) encodes in ([78]) and more generally in ([73) the relative 
positions on the Poincare sphere of the successive polarizer vectors p 



C. Discussion 



The result (|60| ) can be illustrated on the upper sheet H + of an hyperboloid, or, more easily, on its stereographic 

projection known as the Poincare disk. Let us consider in Fig. || the hyperbolic triangle ABC, the three sides a,b, 

and c, of which corresponding to polarizers Pi, P2 and to the resulting polarizer P respectively. In that triangle, 

7i 72 7 

the lengths of sides a,b, and c, opposite to vertices A,B, and C are tanh — , tanh — and tanh — respectively. They 

physically correspond to the quality of the corresponding polarizers: the length of each side approaches the limit 1 
when its corresponding 7 tends to 00. In that triangle (see Fig. |^), a is the angle between the axis of Pi and the axis 
of P2, 4> is the angle between the axis of Pi and the axis of the resulting polarizer P, and 9 is the angle between the 
axis of P and the axis of Pi ■ For the sake of simplicity, we choose the axis of reference along the direction of the first 
polarizers, i.e. u = p 1, so that in Fig. R we have cti = and a 2 = a. 
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FIG. 8: The hyperbolic triangle ABC corresponding to the polarizers Pi, P2 and to the resulting polarizer P in the Poincarc 
disk. 

Now, in order to understand (|60|), let us use Fig. || when the second polarizer P2 is perfect, i.e. b = tanhy = 1. 
In such a case, (|60|) reduces to 



tanh^+e- 
tanh - = 1 



1 + tanh y e la ' 



(81) 



So, using (p9|) , we obtain 



tanh^e^ = e la e~ lS . 
2 



(82) 



This result shows that the polarizer equivalent to using polarizers P\ and P2 successively (and in that order) is a 
perfect polarizer (tanh — = 1) the axis of which is in the direction 4> = a — 5. That result could be surprising, since 
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we could expect to find its axis in the a direction. This comes from properties of hyperbolic triangles. The angle 9 
between P and Pi can be calculated using (pCl): noting that the angle of AC with CB is —a and taking now the 



72 

direction of the arbitrary reference axis for angles in the direction of P2, we get (note that tanh — = 1) 

tanh + tanh e~ la 1 + tanh ^ e~ ia 

tanh le-* e = ^ \ = 1 = 1 (83) 

2 1 + tanh — tanh — e~ ta 1 + tanh — e~ la 

2 2 2 

so that 9 = 0: the axis of P equivalent to using Pi and P2 successively is, as expected, directed in the direction of the 
axis of P2. These two last results, i.e. the angle of P with Pi is <f> = ct — 8 and that of P with P2 is 9 = 0, come from 
the fact that the sum of the three angles of an hyperbolic triangle is not tt, but tt — 8 where 8 is its angular defect. It 
is interesting to underline that this case (P2 perfect) corresponds to the case of aberration in special relativity. 

For the sake of completeness, let us show that 8 is indeed the angular defect, that we note 8', of the hyperbolic 
triangle ABC. The sum of the angles of the hyperbolic triangle ABC gives 

cf> + (tt - a) + 9 = tt - 8' (84) 

so that 

e iS' = e *a e -<* e -tf ( 85 ) 



Using ( A. 16 ) and extracting the angle 9 from the first equality of (p3|), Eq. ( p5|) can be rewritten as 



1 + tanh — tanh — e ia 
e iS ' = ^ 2 (g6) 



1 + tanh -P- tanh — - e 



let 



Comparing ( pq ) with (|5^) gives us 8' = 8. Thus the parameter 8 of the rotation R-^ s (^5|) is the angular defect of 
the hyperbolic triangle build using the parameters of Pi , P2 and P. 



To conclude this section, let us remark that the angle <j) (A. 16) giving the direction of the resulting polarizer P can 
be rewritten as 



^cuclid 



I (87) 



where 



Euclid = arg (^tanh — + tanh — e la 
would have been the angle between the axis of the polarizers Pi and P if the geometry were euclidean. 



VI. WIGNER ANGLE 

The composition law (62), i.e. the addition law ©, is noncommutative |l3| ]. Indeed, the final polarization of the 
light wave passing through the polarizer 1 then through the polarizer 2 is not the same as the final polarization of 
the light wave passing through the polarizer 2 then through the polarizer 1. Such a noncommutativity is obvious in 
the case of perfect polarizers since the final polarization always corresponds to the polarization vector of the second 
polarizer. With two polarizers with respective absorption terms 71 and 72 and with respective vectors ~p 1 and ~p " 2 , 
the two possibilities of composition are 



r ie2 e*® 2 = Tie 4 " 1 8 r 2 e iQ2 (89) 



J4>1<S 
X lffl2 r 

and 

T^e^ 2 ® 1 = T 2 e M2 © Tie lai . (90) 
From the definition of the addition law (|6l]) we easily see that 
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FIG. 9: Representation of the Wigner angle w (194) as a function of a for different value of the product vTTrJ from 0.6 (dot 
line), 0.8 (dot-dot-dash line), 0.95 (dot-dash line)70.99 (dash line), 1 (solid line). 



and consequently 



Tie 



(92) 



However the angles 0i©2 and 02©i are not equal. Thus, by analogy with special relativity (or with optical studies of 
multilayers |p^| ), we introduce the Wigner angle w = </>2©i — ( / ) i©2 which can be seen as a measure of the noncommu- 
tativity of the addition law 0. Using (Elj), ( |89| ) and (|90|), we easily determine the Wigner angle Q (sec also Q and 
@) through 



which gives then 



r 2 e M2 e rie Mi _ 1 + r^e' 

rie iQl ©T 2 e M2 ~ i + r 1 r 2 e- 



2arg(l + r 1 r 2 e J 



rir 2 e 4 



(93) 



(94) 



where a = a 2 — «i is the angle between p 1 and "p^ 2 - 

For perfect polarizers = 1, the light wave is polarized first in one direction ~p \ then in the other direction ~p 2 i 
so the Wigner angle is w = a. For non perfect polarizers the measure of the Wigner angle w can be seen as a measure 
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of the quality of the polarizers. Indeed, if one measures w < a for a given a, we are able to retrieve the product 
and then to state on the quality of the polarizers. Fig. || shows the Wigner angle w as a function of a for different 
value of the product Ti^. 

VII. CONCLUSION 

Textbooks usually only consider perfect polarizers (i.e. polarizers with 7 — ► 00) and do not examine the question 
of knowing how the polarization of a light beam is progressively transformed when propagating inside a polarizer. In 
the present paper, we have addressed these two problems using the notion of Stokes vector and that of the Poincare 
sphere in order to characterize the evolution of a light wave polarization through a polarizer. Whereas the action of a 
polarizer is intuitively thought of to be a simple rotation of the polarization along a geodesic on the Poincare sphere, 
we have shown that it is not always the case. In the general case, more complex trajectories differing from geodesies 
can in fact be found, especially in the case of total reflection based polarizers (see section |EII C ) . How these trajectories 



will affect geometric phases of light wave going through polarizers will be a question addressed in a following paper 
& 

Non perfect polarizers, i.e. here polarizers with finite 7, lead us naturally to define a generalized Malus law and a 
degree of polarization r)(z). 

The action of a polarizer on light wave polarization can be described by Lorentz operators, where the attenuation 
factor 7 is the counterpart of the rapidity $ commonly used in special relativity. The association of two such polarizers, 
involving then two Lorentz boosts, can always be viewed as the result of a Lorentz boost transformation followed 
by a rotation. In other words, the association of two polarizers is equivalent to the association of a polarizer and 
for example a device with optical activity giving an additional rotation of the polarization. We have then derived 



a composition law (61) for polarizers similar to the composition law for noncollincar velocities in special relativity. 
This composition law © appears to be a very natural law to "add" quantities the values of which arc limited to finite 
values. In fact, equation ( pi] ) implies that no matter what are the values of the complex quantities a and 6, such as 
\a\ < 1 and |6| < 1, the modulus of the overall resulting quantity a© b cannot exceed unity. Because it avoids infinity, 
such a generalization of Einstein's composition law of velocities appears to be a natural "addition" law of physical 
quantities in a closed interval. 

The association of N polarizers can also be viewed as the association of a resulting polarizer and an optical device 
giving an additional rotation. Recursive iterations of the composition law (|6l| ) give directly the characteristics of the 
resulting polarizer, i.e. 7 and p , and of its associate rotation, i.e. 5 and q . We have defined the Malus law for the 
association of N polarizers and the corresponding degree of polarization, the latter being easily determined using the 
composition law. 

In addition to the degree of polarization (or equivalently the extinction ratio) we have defined others quantities 



characterizing the quality of the polarizers. We have defined, see section IV, the angle a(z) measuring how far is 



the light wave polarization vector from the polarizer vector p . Also, using the existing isomorphism between special 



relativity and polarization in optics, we have defined, see section VI, the relativistic Wigner angle for the association 
of two polarizers. The Wigner angle is a direct consequence of the noncommutativity of the addition law © ( |6l| ) and 
we have shown that the latter can be used to directly measure the quality of the association of two polarizers. 

APPENDIX 

We aim to determine the absorption term 7, the angle S, and the vectors ~p and ~~q entering the relation (|55|), i.e. 

B^ 2 ,~t2-B-p ll7l = R-q",sB^ n . (A.l) 
We define a = (~Pi,~P2) € [0)""] as the angle between the two polarizer vectors ~p \ and p 2 - Using the definitions 



(W2|) and (U3j) and expanding the exponential operators in (A.l) we obtain the following set of four equations 



72 7i 72 7i $ 7 

cosh — cosh h sinh — sinh — cos a = cos - cosh — (A. 2) 

2 2 2 2 2 2 v ; 



5 7 j, 

sin — sinh — q ■ p = (A.3) 
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sinh — cosh — "oo + sinh — cosh — l?i = sinh — cos — Id — sin — sinh — ~q A ~v 
2 2 2 2 2 2 2 2 



(A.4) 



72.,7i_> A _> <5 , 7 — > 

smn — smn — p 2 A p 1 = sin — cosn — g 



(A.5) 



1. Rotation parameters: <5, g 



As the equation (A. 3) holds for any parameter 7 G IR + and f or a ny angle 5 s [0, 27r[, wc have q ■ p 



Consequently, the two vectors p and q are orthogonal. Equation ( |A5| ) fixes then the unit vector q as 

-» 1 _> _> 
■ p 2 A pi. 



sma 



(A.6) 



The vectors ~pi,~p2 and ~p are coplanar since ~p" is orthogonal to ~q (see Fig. ^). Using (A.6), equation ( A.5 ) can 
be replaced by the following equation 



2 2 

With equations flA.2j) and (|A.7|) the angle 5 is given by 



smn — smn — sin a = sin — cosn — . 



(A.7) 



tan 



7i 72 
r tanh — tanh — sin a 
_ 2 2 

^ 1 + tanh -i tanh cos a 



(A. 



or equivalcntly by 



— = arg ( 1 + tanh — tanh — e 4 
2 V 2 2 



(A.9) 



Thus, from the initial parameters, 71, 72, p 1 and p 2 , equations (A.6) and (A.9) can be used to obtain respectively 
the axis ~q and the angle 5 of the rotation entering (A.l). 



2. Boost parameters: 7, p 



Let us now define the angles ot\ = (~u , ~pi) and «2 = (~u, ~P2) (see Fig. where is an arbitrary reference axis 

in the plane defin ed by p\ and ~pi- In order to calculate the angle <fi = (~u,~p), we project the two members of 
the equation ( |A.4j ) along the direction ~u and along the orthogonal direction ~u A~q. We obtain the following two 
equations 



72 7i .71 72 <5.7 S . 7 

sinh — cosh — cos a? + sinh — cosh — cos a-t = cos — sinh — cos d> — sin — sinh — sin < 
2 2 2 2 22 22 

72 71 .71 72 <5 . 7 .5.7 

sinh — cosh — sin ao + sinh — cosh — sin ai = cos — sinh — sin d> + sin — sinh — cos <t 
2 2 2 2 22^22^ 



(A.10) 
(All) 



The combination of these two equations gives us 

sinh ^ cosh l^e 1 " 2 -f 



sinh — cosh — e l 



2 



Combining now (|A.2| ) and ( A7), we obtain 



(A.12) 



cosh cosh 
2 2 



sinh — ^ sinh ^-e l 
2 2 



coshp*/ 2 . 



(A.13) 
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Dividing the members of equation ( A.12| ) by those of equation ( A.13| ), we find 

tanh e ia2 



tanh - = 
2 



tanh ^ e lai 



1 + tanh tanh ^ e 1 ' 02 "" 1 ' ' 



7i 
2 



(A.14) 



From the initial param eters , 71, 72, 1?! and "p^ 2 > equa tion ( J A.14 ) allows us to calculate the absorption term 7 and 



the direction p entering (|A.l[ ). Indeed, the modulus of (|A.14[) gives us 7 since 



7 

tanh — = 



tanh — + tanh — e i(a2 ~ ai) 
2 2 



1 + tanh — tanh — e l(Q2 - Ql) 
2 2 



(A.15) 



and the direction of the unit vector p is determined from 

7i 

; iami 

a x + arg 



tanh 



' 2 gi(c(2 — «l) 



1 + tanh tanh — e^"^) 



7i 
2 



(A.16) 
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